pub fn new_birthday_probability(n: u32) -> f64 {
    let mut ans = 1f64;

    for i in 0..n {
        let tmp = (365 - i) as f64 / 365f64;
        ans *= tmp;
    }
    1f64 - ans
}
